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Abstract 

00 

^vq The generalised eigenvalues for a pair of TV x TV matrices (Xi,X 2 ) are defined as 

the solutions of the equation det(Xi — XX 2 ) = 0, or equivalently, for X 2 invertible, as 
the eigenvalues of X 2 ~ 1 Xi. We consider Gaussian real matrices Xi,X 2 , for which the 
generalised eigenvalues have the rotational invariance of the half-sphere, or after a frac- 
tional linear transformation, the rotational invariance of the unit disk. In these latter 
variables we calculate the joint eigenvalue probability density function, the probability 
PN.k of finding k real eigenvalues, the densities of real and complex eigenvalues (the 
latter being related to an average over characteristic polynomials) , and give an explicit 
'— 1 Pfaffian formula for the higher correlation functions P(k 1 .k 2 )- A limit theorem for PN.k 

is proved, and the scaled form of p(k 1 .k 2 ) is shown to be identical to the analogous limit 
for the correlations of the eigenvalues of real Gaussian matrices. We show that these 
correlations satisfy sum rules characteristic of the underlying two-component Coulomb 
CO gas. 
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1 Introduction 



The general topic of our study is the statistical properties of the eigenvalues of N x N non- 
symmetric random matrices with real entries. Such matrices will, in general, have both real 
and complex eigenvalues. In the case that the entries of the matrix are independently chosen 
as standard Gaussians — referred to as the real Ginibre ensemble after |20j — a result of 
Edelman et al |13| tells us that the expected number of real eigenvalues is asymptotically 
equal to ^/2N/tt. Numerical evidence presented in the same paper indicates that this 
asymptotic value persists with the standard Gaussian replaced by any distribution of zero 
mean and unit standard deviation. 

Similar properties hold true of the generalised eigenvalues of a pair of N x N random 
matrices (X\,X2) with standard Gaussian entries. The generalised eigenvalues A are spec- 
ified as the solutions of det(Xi — AX2) = 0, or equivalently, as the eigenvalues of X% X±. 
The statistical properties of the A are the specific concern of this paper. Our starting point 
will be to first establish the joint matrix distribution of Y = X^ l X\, which we show in 



Proposition 2.1 to be the matrix Cauchy distribution 

T(Y) _ ^ H n %j;l%i /2) + YY*)->. (1) 

Studies into the statistical properties of the A were initiated in [13] using a different 
logic. It was shown, for example, that the expected number En of real (generalised) 
eigenvalues has the exact evaluation 



^T((N + l)/2) ^ nN 
N T(N/2) tv->oo V 2 1 ; 

To derive ^ , the generalised eigenvalue problem in the case of X\ , X2 having standard 
Gaussian elements was placed in the context of integral geometry. 

First the pair of matrices (XifX?) can be regarded as two vectors in ~R N and the 
corresponding plane spanned by these vectors intersects the sphere S ~ 1 to give a great 
circle. The real generalised eigenvalues correspond to the intersection of this great circle 
with the set A^r of all N x N singular matrices X such that TtXX T = 1 (thus choose 
X = c{X\ — AX2) for suitable c). With X\,X% having standard Gaussian entries, the 
great circle has uniform measure, so the expected number of real eigenvalues is equal to 
the expected number of intersections of A^r with a random great circle. 

Another feature of the random generalised eigenvalue problem studied in [13] is the 
density P(n(A) of real generalised eigenvalues. By writing A = tan# the generalised eigen- 
value equation reads det(cos6>Xi — sinflJ^) = 0. Using the fact that a pair of standard 
Gaussians {x\,X2) is, as a distribution in the plane, invariant under rotation, it was noted 
that (cos#, sin#) must be distributed uniformly on the unit circle, and so 

P(i)(A) = i T f^. (3) 
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The transformation A = tan# is the stereographic projection of the real line on to a 
great circle of the sphere. Prom the famous circle theorem \21\ EJ |37j relating to eigenvalue 
densities of large N x N matrices with general iid entries drawn from any distribution 
with mean zero and fixed a, one might similarly expect uniform asymptotic density on the 
sphere for A~ 1 B where A, B also have iid entries from any distribution with zero mean 
and fixed a — a kind of spherical law^} Certainly, in the case of complex Gaussian entries, 
this is true for every value of N. Indeed, with 

C := -f3B + aA, D := aB + f3A, (4) 

where a,/3 G C such that \a\ 2 + |/3| 2 = 1, one has that (A,B) has the same distribution 
as (C,D), implying that the corresponding joint eigenvalue distribution is invariant under 
rotation of the sphere |26j . In the traditions of random matrix theory — for example 
the circular ensembles of Dyson — the present real Gaussian matrices A _1 B may be said 
to form the real spherical ensemble. In [26] the complex Gaussian matrices A~ 1 B were 
referred to as the (complex) spherical ensemble. 

From ^ and the fact that with A = tan 9, d\/{l+ A 2 ) = d9,we see that as a distribution 
on the sphere, the density of real eigenvalues is uniform. One can also stereographically 



project the complex eigenvalues on to the sphere; it will turn out (eq. (59) below) that in 
the large N limit the total eigenvalue distribution is uniform. Thus, to leading order, the 
concentration of real eigenvalues on a great circle of the sphere does not affect the overall 
eigenvalue distribution. This is analogous to the situation for the real Ginibre ensemble, for 
which the expected number of real eigenvalues is asymptotically y/N/2ir, yet the eigenvalue 
density is to leading order uniform Also, we remark that there is an analogy with the 
random polynomials 

n / \ 1/2 

Pn( Z ) = Y,{) Q P ZP > <*P ~ #[0,1]- ( 5 ) 

p=o 

When stereographically projected onto the sphere there is of order y/~N zeros on a great 
circle corresponding to the real axis [T2], but for N large the density is asymptotically 
uniform on the sphere |31j . 

So in summary, from the pioneering study of the Gaussian real generalised eigenvalue 
problem of [13] , we know the expected number of real eigenvalues and the density of real 
eigenvalues. Beginning with ([!]), this statistical knowledge will be greatly extended. 

An essential ingredient in our analysis is to conformally map from the real line and 
upper half-plane (i.e. the domain of the independent eigenvalues of Y) to the boundary of 
the unit disk by the fractional linear transformation 

z = -. (6) 



This has recently been proved by Bordenave 
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Figure 1: A plot of the eigenvalues from 120 realisations of 100 x 100 matrices with the 
eigenvalues stereographically projected. The great circle of real eigenvalues can be clearly 
seen. 



For A a real eigenvalue, and with e := e te , (|6l) reads 



A = tan? =*(l --?-). (7) 



2 i \ e+ 1 

A key feature is that the corresponding real eigenvalue density is uniform in 9 

= fr> (8) 

as follows from Q. In terms of ^ and Q we show that the explicit form of the joint 
eigenvalue probability density function (jpdf) for k real eigenvalues is 



k (N+k)/2 

Q(Y) = A k<N l[r(e 3 ) J[ rfK)r - A e,w,- , (9) 

j=l s=k+l 1 S| 
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with e = {ei, ...,e fc },w = {wt, wi, W(jV-fc)/2, ™(N-k)/2} and 

(_l)( Af - fc )/2(( Ar -* : )/2-l)/2+(A r -fc)fc/2-fc(fc-l)/4 ^ 



A? 



r(x) 



2 (N(N-l)+k)/2 

xF({N + l)/2) N / 2 T{N/2 + 1) N / 2 , 
x 



\ r(i/2) 2 



^(N-l)/2 









dt 



N/2+1 



1/2 



A e, w, — 
w 



n< 



k (N+k)/2 , 



j<P j=l s=k+l 



W s 
(N+k)/2 



Y[(w b -w a ) (— - — j TT ( — -ti/c), 

ct<£> V S S/ c,d=fc+l V a 7 



which is the content of Proposition 2.3 



We use (|9j) to deduce exact statistical properties. As an example, let PN,k denote the 
probability that there are exactly k real eigenvalues, and set 



TV 



PN,k, 



k=0 



where the asterisk indicates that the sum over k is restricted to values with the same parity 
as N. In Proposition |3.3[ with N even, we show 



z N (0 



(_ 1 )(JV/2)(JV/2-l)/2 



2 N(N-X)/2 
N , N/2-1 



T({N + l)/2) N / 2 T(N/2 + l] 



N/2 



Uwim II 



where 



Oil 



1 r(s/2) 2 

=1 w ' 1=0 

2tt T((A + l)/2) 



(10) 



N-l-Al T(N/2 + 1) ' 
N - 1 — 41 ~T(N + 1) 



N T(2l + l)T(N-2l) r((AT + l)/2) 

^ ; ; A/VT- 



r(Jv/2 + i) 



(ii) 



The product form (10) for should be contrasted with the simplest known forms in 



the case of the real Ginibre ensemble [2 [18], which are determinants of size N/2. 
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It follows from (10) that 



En 



0_ 



z N (0 



2ai 



JV/2-1 



(12) 



where we have used the fact that Zjv(l) = 1- Substituting (11) then reclaims ^ with N 
even. In Corollary 3.9 we give the odd analogue of (12) and with (11) we again reclaim 
(§. The variance in the distribution of the number of real eigenvalues is, by definition, 



given by <j 2 n = ^fc=o k 2 PN,k — E N , which, in terms of the generating function, reads 



2 



+ E N 



E N . 



Substituting (10) and (12), we then obtain the explicit evaluation 

S(N- l/2)r((JV+ l)/2) 2 



2 



2En 



7T- 



T(N)T(N/2Y 



(13) 



(14) 



(2 — t/2)En, which is the same as for the real Ginibre ensemble 



This has large N form aj^ 

m 

Recalling (10) we can read off the explicit value of pN,k- For N, k odd it is a rational 
number, while for N, k even it is of the form Y2i=k/2 c N,2l^ 21 with cn,21 rational (Proposition 

(15) 



3.10). In particular, for N = 2, 



P2fl 



7T 

4' 



and for N = 3, 



P3,i = P3,3 



P2,2 

_ 1 
" 2' 



7T 

4' 



(16) 



As decimals (15) reads p2,o ~ 0.785, ^2,2 ~ 0.214. These values, further approximated to 
0.79 and 0.21 respectively, have been reported in a simulation of the N = 2 case [2H]. The 
latter study was motivated by the question of determining the typical rank of a 2 x 2 x 2 
array (tensor) (see, e.g., [29]). 

The meaning of a tensor in this setting relates to structures A = (a-ijk) G M pxpx2 
represented as the column vector vec A G M 4p2 . As reviewed in [27] , a point of interest is to 
find matrices U = [u x ■ ■ ■ u R ] G W xR , V = [vi ■ ■ ■ v R ] G R pxR and W = [wi ■ ■ ■ w R ] G 
such that 



o2xR 



vec A = ^2 ® V r ® 



r=l 

for i? as small as possible. The positive integer i? is referred to the rank. With both 
( a iji) =: -Xi G M pxp and (0^2) =: -^2 G M pxp random matrices, entries chosen from a 
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continuous distribution, one has that R = p if all the eigenvalues of X 1 X2 are real, and 
R = p + 1 otherwise. 



The simple structure of ( 10 ) also allows for the computation of the large ./V form of the 



distribution of the number k of real eigenvalues. In Proposition 3.5 we prove that it is a 
standard Gaussian in the scaled variable (k — En)/<jn. Further, our method of derivation 
of (10), which involves first computing the exact functional form of the eigenvalue jpdf ([9]), 



allows for the result ^ to be generalised. Thus, in Theorem 4.1 we give a (fci+A^) x (fci+fo) 



Pfaffian formula for the exact (k\, A;2)-point correlation function between k\ real eigenvalues 
and k2 complex eigenvalues. The simplest case beyond ^ is the density of the complex 
eigenvalues in terms of the coordinates Q. With r := \ w\ we find that the complex density 
P(i)( w ) depends only on r, and is given by 

c , , N(N-l)n \ N - 2 /l \ r° dt 
P^ w )= a^+i^ {r +r J Kr- r ) Az^ (l + ^)W ( 17 ) 



The densities (17) and pi) are related by the sum rule for the total number of eigenvalues 



J J d6 p r (1) (0) + 2 J 1 rdr d6 p c (1) (w) = N, 

where the factor of 2 is required since p9^\ only refers to one of each complex conjugate 
pair. 

Section 4 continues with the evaluation of the average of two characteristic polynomi- 
als in terms of elements of the correlation kernel. We also compute the scaled limit of 
the (k\, &2)-point correlation function, and demonstrate that it agrees with the recently 
obtained [9] scaled correlation function for the bulk eigenvalues in the real Ginibre ensem- 
ble. An analogy with a Coulomb gas allows for the formulation of sum rules relating to 
the screening of the effective charge of a fixed number of eigenvalues, and also allows us 
to isolate a certain combination of one- and two-point correlations for which the complex 
moments vanish. 



2 Joint probability density functions 
2.1 Element distribution 

For A, B matrices of size N x N taken from Ginibre's real ensemble so that 

V{A) ■ V(B) = {2^)- N \^ AAT+BBT \ (18) 

we wish to express the probability density function V(Y) of the elements of the matrix 
Y = A~ l B in terms of the elements of Y. This we read off from the calculation of 
V(Y)(dY), where idY) is the wedge product of the independent elements of the matrix 
dY. 
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Proposition 2.1. Let A,B be N x N matrices drawn from the real Ginibre ensemble, and 
let Y = A~ 1 B. The density function for the distribution of the elements ofY is then given 

by®. 

We delay the proof of this result until Appendix [Aj 

2.2 Eigenvalue jpdf and fractional linear transformation 

For a general N x N non-symmetric real matrix, we will have < k < N real eigenvalues, 
where N has the same parity as k. From knowledge of ([!]) we can extract the eigenvalue 
distribution for each allowed k. In this task we are motivated by the work of Hough et 
al [23] (see also and |16j). In particular, we work with the real Schur decomposition 
Y = QRnQ t , where Q is real orthogonal (each column is an eigenvector of Y, with the 
restriction that the entry in the first row is positive) and 



Rn 



Ai 



Rl,k Rl,k+1 

Afc Rk,k+i 

Zk+l 



R 



l.m 



Rk,m 
Rk+l,m 



m = (N + k)/2. 



(19) 



Rn is block upper triangular: on the diagonal we have the real eigenvalues Aj and the 2x2 
blocks 



.X j 



bj,Cj > 0, 



which relates to the complex eigenvalues Xj ± iyj , yj = WbjCj. Note that the dimension of 
Rij depends on its position in R^: 

• 1 x 1 for i, j < k, 

• 1 x 2 for i < k, j > k, 

• 2 x 2 for i,j > k. 

For this decomposition to be 1-to-l we need the eigenvalues to be ordered and we choose 

Ai < • • • < Afe and xt+i < • • • < x m . (20) 

Since we are looking to change variables from the elements of Y to the eigenvalues of 
Y, as implied by the real Schur decomposition, before proceeding we first need knowledge 



S 



of the corresponding Jacobian. Prom [TT] we know that 
{dY) = 2( JV - fc )/ 2 niA( J R PP )-A( J R ji )| 

k (N+k)/2 

x(dR N )(Q T dQ) Y[d\ s \b t - c t \ dxidbidci, 

s=l l=k+l 

where \(Ru) = A; for / < k, while \(Ru) = x\ ± iyi for I > k, and Rn is the strictly upper 
triangular part of Rjy. Our interest is only in the eigenvalue dependent portion so we can 
immediately dispense with the dependence on Q by integrating out (Q T dQ) according to 



/ 



N 



(see e.g. [33, Theorem 2.1.15], with the modification of omitting the factor of 2^ since we 
have specified the columns of Q to have positive first entry). 

For our goal of computing the eigenvalue jpdf the objective now is to integrate over all 
the Ri,j] this procedure is similar to that in [26] where that author was concerned with the 
analogous complex ensemble. The details in the present setting are contained in Appendix 
[B| According to this working ([I]) has been reduced to the distribution of {Aj, Xj, bi, Cj}, 

n -W-k)/4 r u N + i)/ 2 )^ 2 r(iV/2 + l) N ' 2 ( r (( iV + 1 )/2) \ m rr 1 _L 

V r(JV/2 + i) J y r(j/2) 2 

(N+k)/2 k (N+k)/2 

x s n +l det ( i + zszi)^ n (1 + A? )(- +1 )/ 2 2(JV " &)/2 n +i i&<-«i 

X j^J |-Rpp — -Rjjr'l) (21) 

where use has been made of the simplification 

"tt 1 (k -s)/2 r((iv + i)/2) JV _ 2a _ 2 r((iv + i)/2) 1W2 + 1) 

11 r((7v + fc- s + i)/2) JUL r(JV-a-i/2) r(N-s) 



n (k+N 2 -2N)/4 



T((N + l)/2)^T(N/2 + 1)^ ( ^0^ ) ^ n 



li r((Jv + i + «)/2) 

By writing e s = 1 + x 2 + y 2 s and <5 S = 6 S — c s we see that det(l2 + z s z^) = e 2 + 5 2 . Also, 
from [TT], we know that 

dx s db s dc s = — ^ s :dx s dy s dd s , 
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although we require — oo < 5 < oo rather than < 5 < oo as claimed by Edelman. So now 
we integrate over 5 



»<5=oo 



det(l + Zszf)"/^ 1 



dxedbcdcc 



4y s 
4y s 



<5=oo 

<5=0 
t=oo 



SdS 



( e 2 + <5 2 )^/ 2 + VV + <* 2 
dt 



dx s dy s 



t=2ys (^-4^+^2)^/2+1 



dx s dy s . (22) 



Substituting (22) in (21) as appropriate gives the reduced jpdf, but ([22]) as written appears 
intractable for further analysis. On the other hand it follows from the analysis relating to 
Q with a, P 6 E that when projected on to the sphere the eigenvalue density is unchanged 
by rotation in the XZ plane, where X, Y, Z are the co-ordinates after stereographic pro- 
jection. This suggests that simplifications can be achieved by an appropriate mapping of 
the half-plane that contains the rotational symmetry of the half sphere. 

We therefore introduce the fractional linear transformation (|6j) mapping the upper half- 
plane to the interior of the unit disk, and ([7]) mapping the real line to the unit circle. In 
particular, the complicated dependence on x s ,y s in (22) is now unravelled. 



Lemma 2.2. Let e, 



4y, 



t=oo 



1 + x 2 s + y 2 . With the substitutions and ^ 
dt 



2ys (e 2 _ 4 y2 + t 2)AT/ 



2+1 



dx s dy s 



?/>.. 



2 2N - 2 \w s \ N+1 



dt 



N/2+1 



dUedVe 



(23) 



Proof: Noting that 



1 



w. 



Us 



dx s dy s = 
reduces the given expression to 



16KI 
\l + w s \ 4 
dz- 2 



ducdv q 



|l + ^sl 4 



dUedVe 



16 



1 



w. 



dt 



1 + w s \* |1 + w s \ 2 J 2 (i-\ Ws \2)/\i+ Ws \z (16(M 2 /|1 + u; s | 4 ) + t 2 )*/™ 



duedve. 



The RHS of (23) results from this after the change of variables t *-+ 4\w s \t/\l + w s 



□ 
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For the remaining terms in (21), the substitutions ^ and Q give 

k 



3<P 



s — 1 1 1 j<p 



k {N+k)/2 

n n \ x j - z *\\ x o 

j = l s=k+l 



(_ 1 jk(N-k)/2 2 (N-k)k TT^(iV-fc)/2 
J'=l 



e 7 - + 1 



(N-k) 



(N+k)/2 
s=fc+l 



1 



W s + 1 



2fe fe (N+k)/2 

n n ( 

j = l s=fe+l 



to. - e, 



n 

fc+l<a<fc<(AT+fe)/2 



(N+k)/2 [N+k)/2 

z a -z b \\z a -z b \ J] k-Zdl = (-2) 2 (^^^) JJ (l-l^l 2 )- 1 



c,d=k+l 
c^d 



(N+k)/2 

X JJ (w] 

s=k+l 



JV-fc-1 



1 



IU S + 1 



j=k+l 



HN-k-2) (1 ^ j 

Wfe - t£?o) — Wc 

\w b w a J /J- Vw d 



a<f) 



An essential feature is that, apart from the creation of some one-body terms that depend 
only on the radius, the product of difference structure is conserved by the substitutions. 
Combining all this, and noting 



N-l 



N-l 



1 



cos(0j/2) 



N-l 



we have the explicit form of the eigenvalue jpdf in the variables ^ and Q. 

Proposition 2.3. Let A,B be N x N matrices drawn from the real Ginibre ensemble, and 
let Y = A^ 1 B. In the variables and |?|) the eigenvalue jpdf ofY, conditioned to have 
k real eigenvalues (k being of the same parity as N ) is 



3 Generalised partition function 

3.1 N even 

Having established the eigenvalue jpdf, we wish to find the generalised partition function 
from which we can calculate probabilities and correlations. For {ei, • • •, et\ corresponding 
to real eigenvalues and {wk+i, • • ',iV(N+k)/2> Wk+i-, ■ ■ W(N+k)/2} corresponding to complex 
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conjugate pairs, define the generalised partition function by 

r0a r 6 3 r 27T k 



z HN . k)/2 [uM = {{N _ k)/2)] J o ^\ 6x de 2 ...J &k je k ^u( ei ) 

, , (N+k)/2 

x / dw k+1 --- / dw {N+k)/2 ]"] v(w s )Q(Y), (24) 
Jn Jn 



where the factor of l/((-/V — k)/2)\ comes from relaxing the ordering constraint on the 
complex eigenvalues, and Q is the unit disk. Note the ordering of the angles corresponding 
to the real eigenvalues, in accordance with the ordering of {Aj} in (20). 

It is at this point that parity considerations become important. For the time being, we 
will assume that N (and consequently k) is even. 

Proposition 3.1. The generalised partition function in the case where N is even is 

(--\)(N/2){N/2-l)/2 

Z k , {N - k)/2 [u,v] = { J 2jV(jV _ 1)/2 T((N + l)/2) N / 2 T(N/2 + l) N / 2 

N 1 

x IIf^727 [cfc]pf [C 2 «;,*[M + /^M , (25) 

with [C k ] denoting the coefficient of C k , and where 



otj,k[u,v] = -- / d0 1 u{e 1 )T{ei) I d0 2 u(e2)T(e2)qj-i(ei)qk-i{e 2 )sgn(02 - 9i), 
1 Jo Jo 

Pj,k[ u > v ] = J dw v(w)t(w)t i fq j _i(w)q k -i - q k - 1 (w)q j - 1 (^j) > ( 26 ) 

and 

q 2j {x) =p 2j (x), q 2 j+i{x) = PN-i~2j(x), (27) 

with Pi{x) an arbitrary monic polynomial of degree i. Equivalently, with 

N 

Z N [u,v] := ^ z k,(N-k)/2[u,v], (28) 



fe=0 
k even 



we have 



(--\\(N/2)(N/2-l)/2 

Zn[u,v] = { } 2N(N _ l)/2 r((AT+ l)/2)^ 2 r(iV/2 + 1)^ 

N I r n 

x II r(s/2) 2 Pf r 3 "''^' ^ + ' ^ 
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Proof: With Pi(x) an arbitrary monic polynomial of degree i, the Vandermonde product 
in Q(Y) can be written 



A I e,w, 



det 



[Pl-l( w s)} s=k +l^...,(N+k)/2 

\pi-i(l/w s )} s=k+1 ^ t{N+k)/2 



(_l)(W-fe)/2((A r -fe)/2-l)/2 det 



1=0,... ,N 

\pi-i(ej)]j=i,...,k 
pi-i(l/w s ) 



s=k+l,...,(N+k)/2 



1=1, ...,N 



where, for the second equality, we have interlaced the rows corresponding to complex 
conjugate pairs; this will be convenient later. 

Next we apply the method of integration over alternate variables |32l HH] 



'k,(N-k)/2[U, v \ 



= (-l)( JV - fe )/ 2 (( JV - fe )/ 2 - 1 )/ 2 - 



A 



k,N 



2tt r 2ir 



d9 2 / d9 4 



x / dwk+i 
Jn 



(fc/2)!((iV-fc)/2)!j 'J 

(N+k)/2 



2tt 



dO k 



dw 



(N+k)/2 



n t~12 t ( w ^ t [— 



B=k+l 



xdet 



u(e)r(e) Pl ^(e)d0 
u(d 2 j)r(e2j)pi-i(e 2 j) 



] = l,...k/2 



Pl-l(Ws) 

pi-i(l/w a ) 



=fe+l,...,(iV+fc)/2 



(=l,...,jV 



Re-order columns in the determinant according to 

P0,PN-l,P2,PN-3, • • -,PN-2,Pl, (30) 

introducing a factor of (— l)( JV / 2 )(-^/ 2-1 )/ 2 . For labeling purposes we define 
q2j(x) = P2j(x), q2j+i{x) = PN-l-2j(x). 

Expanding the determinant according to its definition as a signed sum over permutations, 
then performing the remaining integrations gives 



Z k , (N - m [u,v] = ( _l)(^)/ 2 ((^)/2-D/2 



Ak,N 



(fc/2)!((JV-fc)/2)! 

fc/2 N/2 

x ^2 € ^ II a P{2l-l),P(2l) [ J &P(2Z-1),P(2Z)> 
PeSjv Z=l Z=fc/2+l 
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where 



2vr i-6i 

a jk = I d0 1 «(0 1 )r(ei)oj_i(ei) / d02u(O 2 )T(e 2 )qk-i{e2) 



o 



bjk = [ dw v(w)t{w)t ( — ) qj-i(w)q k -i [ 4 



If we now impose the restriction P(2l) > P(2l — 1), (/ = !,..., iV/2) this can be rewritten 



as 



Z HN - m M = (-l)^- fe )/ 2 ^- fc )/ 2 - 1 )/ 2 (2z) fe / 2 A fc(JV 

fc/2 AT/2 

X X] £ ^ J\ a P(2l-l),P(2l) II Pp{21-1),P{21)i ( 31 ) 
-PeSjv 1=1 l=k/2+l 

P(2Q>P(2i-l) 



with ajfi,Pj,k given by (26). But for general C = [cjfe]j ) fe=i ) ... i Ar, c^-fc = —Ckj and with iV 
even, 

N/2 

PfC:= ^ e(P)nc P(2 ,_ 1))P(20 (32) 

PeSjv z=i 

P(2i)>P(2i-l) 



allowing the sum over permutations in (31) to be written in terms of a Pfaffian, and (25) 



follows. □ 

With u = v = 1, -^fc,(Af-fc)/2[' u 5 v ] is then the probability of finding k real eigenvalues and 
(N — k)/2 complex eigenvalues, and it can be calculated using the result of Proposition 



3.1 The real Ginibre ensemble permits an analogous formula, which by choosing the 
arbitrary degree i polynomials Piix) to be the monomial x l , further simplifies to involve 
a determinant of size N/2 x N/2. For the problem at hand we can do even better, by 
explicitly constructing polynomials qi{x) which for general £ skew-diagonalise the matrix 



in (25), with the 2x2 diagonal blocks of the form 

. « „. . „n 11 _i_ «„..,„.. „n n 

(33) 



C 2 a2i+i,2i+ 2 [l, 1] + fe + l,2i +2 [l, 1] 

-(C 2 a2i+l,2i+ 2 [l, 1] + /8»fl,2t+2[l, 1]) 



and all other entries 0. 

The desired so called skew-orthogonal polynomials turn out to be quite simple, and it 
was knowledge of these polynomials that motivated the definition of the {qi(x)} in terms 
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of the {pi(x)} in (27). First we define the inner-products 

f2n 



(Pj,Pl) 







2?r 



d0ir(ei) / d8 2 T(e2)pj(e 1 )pi(e 2 )sgn(62- 9i) 



o 



Oij+l,l+l[l, 1], 



(Pj,Pl)c-= [ dw t(w)t ( 4 ) 



Pj(w)pi 



w 



pi(w)pj 



(34) 



(35) 



and we see that the polynomials that skew-orthogonalise (34) and (35) are exactly the 
polynomials that yield the block-diagonal matrix (33). That is, we look for polynomials 
that simultaneously satisfy the conditions 



(P2j,P2l)r = (P2j+l,P2l+l)r = 0, (p2j,P2l+l)r = ~ (P21+1, P2j) r = 5jJ «/, 
{P2j,P2l)c = (p2j+l,P2l+l)c = 0, {p2j,P2l+l)c = ~ (j>2l+l, P2j) c = Pi, 



(36) 



where oti and Pi are given by (11). 



Proposition 3.2. The skew-orthogonality conditions (36) are satisfied by the polynomials 
Pj(x) = xi and thus, according to (27), 



Q2j{x) 



r 2 i 



Q2j+l 



X = X 



N-l-2j 



(37) 



Further, with ctj,Pj as in (11) 

"2j+l,2j+2[l, 1] = «j, /32j+l,2j+2[l, 1] = Pj, 

and consequently 



Pf [C 2 a^[l,l]+/3 jV [l,l]] 7Z=1 



2 N 



N/2-1 

n (c%+a). 

1=0 



(38) 



Proof: The skew-symmetry property a^[l,l] = —ay [1,1], ^[1,1] = — /3y[l,l] can 
be checked by observation, so to establish the result we must show that both otj t i\\, 1] and 
PjlW, 1] are non-zero only for j = It + 1, 1 = 2t + 2 in which case they have the evaluations 
stated. 

From (26), we have 



aj+i,i+i[l, 1] = c 



2?r 



2t 



d g ie ^U-(N-i)/2) / d02 e ie 2 (i-( N -i)/2 



2t 



d0l e ie,Q-{N -i)/2) / d02e ie 2 (i-(N-i)/2 } 
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where c is a constant factor and j = 2j or j = N — 1 — 2j for j even or odd respectively. 
Performing the inner integrals over 62, using the fact that j 7^ (N — l)/2 since j is an 
integer and N is even, we find 



oy +W [l, 1] = c . I - dBi e*<&-™\ (39) 
21 — TV + 1 jo 

which is non-zero only in the case that / = N — 1 — j, or j = 2t + 1,1 = 2t + 2 (for 
1] positive). The evaluation of (39) in this case is straightforward. To obtain 
the conditions on s and t where /3 S) t[l,l] / we repeat the procedure used above by 
writing out t(w),t(w^ 1 ), q2s,Q2t+i- The fact that s = 2j + 1 and t = 2j + 2 for a non-zero 
evaluation is then immediate. 

It remains to show that /32j+i,2j+2[l) 1] = Pj, which requires knowledge of a non- 
standard form of the beta integral. After converting to polar co-ordinates, setting c := \ w\ 2 
and integrating by parts, one obtains 

2tt 3 / 2 r((iV + l)/2) 2*+^ f 1 c 2 i + c N -^- x 
/%+i,2i+2 - N _ 1 _ 4j r(JV/2 + l) + N-l-4jJ (l + c)^ 1 



According to [22], Equation 3.216 (1)], for general a, b such that Re 6 > 0, Re (a — b) > 0, 

r(6)r(a-6) 



f (t^ 1 + t a - b - 1 )(l + t)- a dt 

Jo 



r(o) ' 

and thus, with b = y, a — b = x, is a non-standard form of the beta integral 

The stated formula for 02s+i,2t+2 now follows. □ 

This simplification allows us to give a simple product form for the generating function of 
the probabilities {pN,k}, and to proceed to specify statistical properties of the corresponding 
distribution. 

Proposition 3.3. Let N be even. The generating function for PN,k> the probability of 
finding k real eigenvalues, as specified by 

N N/2 

Zn(0 ■= *W = J2^ kZ ^,(N-2k)/2[iA], (40) 

fc=0 A:=0 

k even 



has the evaluation (10). 



Using (40) we calculate the expected number of real eigenvalues and the variance. 
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Corollary 3.4. With N even the expected number of real eigenvalues is 

N/2-1 



d 



\ - 2ai 



(41) 



which has evaluation given by p|). The variance in the number of real eigenvalues is 



2 



N/2-1 



Z N (0 =i + E N - E 2 N = 2E N - A 



z=o 



which has evaluation \ 



Proof: The second equalities follow from the first and (10), while to obtain ([2]) and (14) 
use has been made of the summations 



N/2-1 

E 

3=0 



N - 1 
2j 



-.7V-2 



N/2-1 

E 

3=0 



N- 1\2 _ 2 2N - 3 T(N - 1/2) 
2j ) ^rT(N) ' 



□ 



As noted in the Introduction, the result (41) was first derived by Edelman et al. [13] 
using ideas from integral geometry. A corollary, also noted in |13j , is that for N — > oo 

1 



En 



ttN 
~2~ 



1 



1 

4iV + 32A 2 



+ 



128A 3 



21 /J_ 
20487V 4 + ViV 5 



It was also remarked in the Introduction that to leading order ( 14 ) implies the variance is 
related to the mean by a 2 N ~ (2 — t/2)En, which coincident ally (?) is the same asymptotic 
relation as found in the case of the real Ginibre ensemble |17j . 

The explicit form of the generating function given in Proposition 3.3 allows for the 
computation of the large N limiting form of the probability density of the scaled number 
of real eigenvalues. 



Proposition 3.5. Let a 2 N and En be as in Corollary 3.4 
part. We have 



and let [■] denote the integer 



lim sup a;6( _ 0O oo) 



a NPN,[a N x+E N ] 



-x 2 /2 



2vr 



0. 



Proof: For a given n, let {p n (k)}k=o,i,...,n be a sequence such that 

n 

Pn{x) = J2 Pn {k)x k 



k=0 
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N, k 




40 k 



Figure 2: A plot of P3oo,fc> that is, the probability of finding k real eigenvalues from a 300 x 
300 real matrix Y = A~ 1 B, where A, B are real matrices with iid Gaussian elements. These 



values were calculated using (10). The solid line is the Gaussian curve from Proposition 



3.5 with a normalising factor of 2 since N and k must be of the same parity. 



has the properties that the zeros of Pn(x) are all on the real axis, and -FW(l) = 1. Let 



Mr 



k=0 



k=0 



and suppose a n — > oo as n — > oo. A local limit theorem due to Bender |7] gives 



\(7 n p n ([a n x + fJL n }) 



-x 2 /2 



2vr 



(42) 



Application of this general theorem to Zjv(£), with £ 2 = x, gives the stated result. □ 



In Figure 3.1 we have used ( 10 ) to calculate the value of k for iV = 300, k = 2, 4, 38 



and overlaid it with the Gaussian curve given by Proposition 3.5, the agreement is quite 
clear. 



3.2 N odd 

Here we will deal with the case of odd. The parity of N has no bearing on ([9]), nor on 
(24). However, the proof of Proposition 3.1 uses the method of integration over alternate 



variables, which crucially depends on the evenness of N as it pairs the real eigenvalues. In 
the case of N odd, one of the real eigenvalues cannot be paired; it is the treatment of this 
eigenvalue that distinguishes the even and odd cases. 
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Proposition 3.6. With otjj [u, v ] and f3jj[u,v] as in (26), the generalised partition function 
in the case of N odd is 



Z k , {N - k)/2 [u,v] = { ) 2N{N _ 1)/2 T{{N+l)/2) N ^T{N/2 + l) N ^ 



where, with e := e ld (recall ^jj), 



■m u \\j=\,...,ii 



vi[u] :-- 



2tt 



and 



qij =P2j, 

Q2j+1 = PN-l-2j, 



12 j — P2j+1, 
Q2j+1 = PjV-l-(2j+l) 

QN-1 = P(JV-l)/2- 



u{e)r(e) qi ^(e)de, 

0<2j<(N- l)/2, 

(N - l)/2 <2j<N-l, 



The analogous definition to (28) for Zn[u,v] where N is odd is 

N 



Z N [u,v] := ^ Z k,N-k/2[u,v], 



fc=i 

k odd 



and 



Z N [u,v] 



_ 1 )(7V-l)/4((7V-l)/2)-l) 



2 iV(JV-l)/2 



T({N + \)/2) N ' 2 V{N/2 + l) N ' 2 



N 

n 



1 



rPf 



[ai d [u,v] + Pi,j[u,v]] ij=1 N [vi[u]] i=1 



N 



-"iMl^i,..., 



N 







(43) 



(44) 



(45) 



(46) 



Proof: Similar to the even case, again with arbitrary monic polynomials {pi(x)}, we 
write the Vandermonde product of Q(Y) as 



A ( e,w, — ] = dct 



[pi-i{ej)} j= i,...,k-i 

\Pl-1.( w s)]s=k+l,...,(N+k)/2 

_ [Pi-i(iM)] s=fc+li ^ !(JV+fc)/2 



^_ 1 ^(W-fe)/2((JV-fe)/2-l)/2 det 



l=0,...,N 

[p/-i(ej)]j=i,...,fc 

Pl-i(w s ) 
pi-i(l/w s 



[Pz-i(efe)] 



s=fe+l,...,(Af+fe)/2 



(47) 



1=1,. ...N 
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where we have moved the row corresponding to the kth. real eigenvalue to the bottom of 
the matrix. This always involves an even number of transpositions so no overall factor is 
required. It is more convenient converting this matrix to Pfaffian form than the equivalent 
matrix where the kth row is not moved. This row corresponds to the single unpaired real 
eigenvalue that must exist in any odd-sized real matrix, a fact which is guaranteed by N 
and k being of the same parity. 



Now we substitute (47) in (24) and apply integration over alternate variables, as in 
Proposition 3.1 to find 



^k,(N-k)/2[u,v] = (-1) 



(Af-fe)/2((JV-fe)/2-l)/2 



.4 



k,N 



< / d0 2 / dO A • • • / d6 k _ 1 / dw k+ i ■■ dw {N+k y 2 TT v(w„) ^t(w s )t ( — 

Jo Jo Jo Jn Jn „_ii_i \ w s\ \ W ' 



((k-l)/2)\((N-k)/2)! 

(N+k)/2 



s=k+l 



x det 



J" 2 ' u(8)T(e)pi-i(e)M 
u(e 2 j)T(e 2 j)pi-i{e2j) 



i=l,...(fe-l)/2 



Pl-l(w s ) 

Pl-x{l/w s ) 



s=fe+l,...,(JV+fc)/2 



l=l,...,N 



We need to reorder the columns of the determinant in a similar way to that of (30) 



although with the key difference of shifting the middle column to the end. As in the even 
case, this is to assist in the use of skew-orthogonal polynomials. The re-ordering becomes 



where 



P0,PN-l,P2,PN-3, ■ ■ ■>P(jV-l)/2-e 1 ,2>P(Ar-l)/2+ei, 2 > 
P(W-l)/2+e 2 ,i>P(JV-l)/2-e 2 ,i) " " 'PN-4, P3, PN-2, Pi, P(JV-l)/2i 



(48) 



e l,2 



£2,1 



1 for (JV- l)/2 even, 

2 for (N - l)/2 odd, 

2 for (JV- l)/2 even, 

1 for (TV- l)/2 odd. 



This introduces a factor of (_i)(^-i)/2+(A f -i)/2((7V-i)/2-i)/2 > AlsQ5 for N odd) the factors 
of (—1) in j4fe jn can be re-written by noting 



(-1) 



(JV-fc)jfe/2-jfc(fc-l)/4 



(-1) 



(JV-l)/2-(fc-l)/4 
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This gives us an overall factor of 

/_^(Ar-fe)/2((Ar-fc)/2-l)/2 x /_ 1 ^(JV-l)/2+(iV-l)/2((7V-l)/2-l)/2 x ^ N 



(_l)(A r -l)/4((W-l)/2)-l)-(fc-l)/4 
2(N(N-l)+k)/2 " 



TV 



((N + 1) /2) N ' 2 T{N/2 + 1) N / 2 l\ 



Now we again expand the determinant as a signed sum over permutations and impose the 
restriction P(2l) > P(2l — 1). This gives us the odd analogue of (31) 



N 



C_i\(AT-l)/4((iV-l)/2)-l) - i 

z k , {N - k)/2 [u,v] = ( > 2N{N _ 1)/2 r((iv + i)/2) w / 2 r(iv/2 + i) w / 2 n 



ir( s /2) 2 



(fc-l)/2 (N-l)/2 

^ e(P) Vp(N),N+l Yl Ctp(2l~l),P(2l) Yi Pp(2l-l),P(2l) 
1=1 Z=(fc+l)/2 



pes N 

P(2I)> P(2I-1) 



where i/p(7\r)[u] := ^p(Ar) ; Ar+i is given by (44). Using the Pfafnan definition (32), (43) now 
follows. □ 



To simplify the calculation of the Pfafnan for the even case in (25) we used the skew- 



orthogonal polynomials of Proposition 3.2 We would like to find equivalent polynomials 



for the odd case, that is polynomials that will reduce the Pfaffian in (43) to the block 
diagonal form 



B 1 

'-. 

••• 



B 



(JV-l)/2 










v N 

-V N 



(49) 



where the Bi are the 2x2 blocks given by (33). However, the best we can do here is to 
obtain the structure 



B L 


-Gx,n 








B 



(JV-1J/2 J 

-Gn,n 










Gn,n _ 



v N 
-i/jv 



(50) 



where G S) n '■= C a s,iv[l>l] + /? Si jv[1> !]• The structure (50) will be sufficient for our ends 



since the Pfaffians of ( |49| ) and (50) are equal, which can be seen by applying the Laplace 
expansion for Pfaffians. 
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By comparing (50) with (33) we see that we are looking to skew-orthogonalise the 
same inner products (34) and (35) as in the even case, however those inner products are 
dependent on N. Also, the column reordering (48) means while these polynomials are still 
monomials, compared to the even case the labeling is more complicated, since there was the 
additional movement of the middle column to the end. The first half of the polynomials are 
the same as the even case, while the second half are modified by j — > j + 1/2. The middle 
polynomial must be singled out for special treatment. For these reasons, the specification 
of the skew-orthogonal polynomials for TV odd is different (and more complicated) from 
that for N even. 

Proposition 3.7. The matrix 

[CV;[M] + /MM]] ii=1 N ^[i]U,...,iv 
[-"iMW,...,* 



evaluates to the modified block diagonal form of (50) using the polynomials Pj(x) 
(j (N — \)/2), and thus according to (45) 



Q2j{x) = 
<72j+iO) 

Q2j(x) = 



X 



,2j+l 



N-l-2j 



0<2j<(N- l)/2, 



(N - l)/2 < 2j < N - 1, 



provided the degree (N — l)/2 polynomial qN-i(x) is chosen as 

(AT-l)/2-l T (JV-l)/a\ 



x 



(N-l)/2 



E 

3=0 



(l2j+l,- 



12 3 + 1,X 



(AT-l)/2\ 



d2j+l,2j+2[lj 1] + /?2j + 1,2j+2[1i 1] 
(q2 3 ,x( N -^) r + (q 2j ,x( N -^) c 



Qf2j+l,2j+2[l, 1] + P2j+l,2j+2[l, 1] 

With these polynomials 



#y+i,2i +2 [l,l]=&, 



Q2j+l( 



< 2j < (N - l)/2, 

(N - l)/2 < 2j < N - 1, 
a S)A r[l, 1] + /3 S)A r[l, 1] = a N> .[l, 1] + ^,.[1, 1] = 0, 



a2j+i,2i+2[l, 1] — a j+1 /2, 

^2j+l,2j+2[l, 1] = /3j+l/2) 



s < N, 



v N ■= v N [l] 



7T 



r((iv+i)/2) 
r(Ar/2+i) ' 
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where ctj,/3j are as in (11) and 



2tt r((JV + l)/2) 
N-S-Aj r(JV/2 + l) ' 



ft 



2aAt 



^•+1/2 - 

So in particular 



4j 



)N r(2j + 2)T(N-2j-l) 
T(N + 1) 



T((N + l)/2) 
T{N/2 + 1) 



Pf 



[^[ 1 ]]i=l,...,JV 





[(HilMRftjlM]]^ 

\- V A l \\j=l,...,N 
r(JV-l)/41-l (7V-l)/2-l 

=vn n (c 2 ««+a) n (c 2 %i/2+a+i/2), 

(=0 i=f(JV— 1)/4"1 

where \x\ is the ceiling function on x. 



(51) 



Proof: For < 2j < (N — l)/2 we have the result by Proposition 3.2 and replacing 
j \— > j + 1/2 we have result for (N — l)/2 < 2j < N — 1. By the construction of qN-i, we 
see that a S) jv[l, 1] + /3s,at[1> 1] = for 1 < s < N. All that remains is to show that the row 
and column of ui obey the skew-orthogonality condition. 

Writing out the z//[l] in full, the fact that it is non-zero only for Z = iV is clear, that is, 
only when I = N does the angular dependence cancel from the integral. In which case the 
evaluation is straightforward. 

Recall from above that the Pfaffian of the modified block diagonal structure (50) is 
equal to that of the Pfaffian of (49) and so we have the evaluation in (51). □ 



The odd analogue of ( 10 ) can now be given. 
Proposition 3.8. In the case of N odd, the generating function for Pn,u is 

N 



Z N (0 := £ tpN, 



k odd 



(JV-l)/2 

£ 2m ^2«+l,(Ar-2«-l)/2[l> 1] 

1=0 



and has evaluation 

(_!)(iV-l)/4((JV-l)/2)-l) 



N 



z N (0 



2 N(N-l)/2 

r(jv-i)/4i-i 



T((N + l)/2) N ' 2 T{N/2 + l) N ' 2 Y\ 



(jV-l)/2-l 

X l\ (£% + A) ] [ (£Vh/2+A+1/2)- 

1=0 Z=[(AT-l)/4] 



(52) 
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From Proposition 3.8 we can calculate the expected number of real eigenvalues in the 



case of N odd, which we know from [13] is given by ^ independent of the parity of N. 
Similarly, we can check that the formula (14) for the variance also holds independent of 
the parity of N. 

Corollary 3.9. For N odd, the expected number of real eigenvalues ofY can be written 

|-(.iV-l)/41-l (JV-1J/2-1 

e n = i+ v _^- + v 2ai+ y , 

which has evaluation M). 

The variance for IS odd is 



„2 



4 a 2 , 

^ u Z+l/2 



r(JV-l)/41-l 4q _ 2 (JV-l)/2-l 

2(En " 1} " § + ^ r( S )/41 W 1/a + a + i/ 2 ) 2 ' 



which has evaluation (14)- 



Proof: The formulae in terms of {ai,j3{\ follow from Proposition 3.8 and for the ex- 
pressions for En, o~ 2 n in terms of Zn(£) recall (12) and (13). For the summations we use 
the identity 

r(JV-l)/41-l , N _ (iV-l)/2-l p (JV-l)/2-l _ ^ p 



Z=0 V 7 i=[(JV-l)/4] V 7 Z=0 

for integer p and for both (N - l)/4 € Z and (JV - l)/4 € Z + 1/2. 



□ 



The values of Pjva for N = 2, 7, calculated using Propositions 3.3 and 3.8 are listed 



in Table [TJ along with the results of a simulation of 100,000 matrices. A remarkable fact 
can be immediately seen in the table: the probabilities for even iV are polynomials in tt 
of degree N/2, while for odd N they are rational numbers. The key difference is that 
(N + l)/2 and N/2 + 1 alternate as integers and half integers, depending on whether N is 
even or odd. These values introduce factors of y^r through the gamma functions. 



Proposition 3.10. Let PN,k be the probability of finding k real eigenvalues in a matrix 
Y = A- X B, where A, B are Gaussian real. Then for N even, pN,k is a polynomial in tt of 
degree N/2. For N odd, PN,k is a rational number. 

Proof: For iV even ai and the second term in /3/ (with u = v = 1) both yield factors of 
7r 3 / 2 . The pre-factor in (25) yields 7r _iV//4 . Combining these two we find the highest power 



of tt is N/2. Noting that the first term in /?/ has a factor of ir 1 / 2 and expanding the product 



in (25) gives lower order terms in tt. 



For the odd case, the pre-factor in (52 ) gives tt 



-N/4-1/2 



. Then by noting that (^ai + fii) 
and (£ 2 a;+i/2 + A+1/2) both give factors of 7T 1 / 2 and vn gives 7r 3 / 4 we see that the end 
result is a rational number. □ 
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Exact Vat b 


Decimal pjq & 


Simulated pjy ^ 


P2,2 




0.785398 


0.78691 


P2,0 


1 - |tt 


0.214602 


0.21309 


P3,3 


i 

2 


0.5 


0.50051 


P3,l 


1 
2 


0.5 


0.49949 


P4,4 


1024 


0.260234 


0.25705 


P4,2 


3^ _ 27 2 


0.65763 


0.66053 


P4,0 


1 _ 3 + _27_ 2 


0.0821365 


0.08242 


P5,5 


1 
y 


0.111111 


0.11167 


P5,3 


11 
18 


0.611111 


0.60969 


P5,l 


5 
18 


0.277778 


0.27864 


P6,6 
P6,4 
P6,2 
P6,0 


84375 „3 
67108864 /l 
14625 2 253125 _3 
262144 " 67108864 " 
15 14625 _2 | 253125 3 
32" 131072'' 1 67108864" 
1 15 , 14625 _2 84375 _3 
x 32" ^ 262144'' 67108864" 


0.0389837 
0.433673 
0.488323 
0.0390194 


0.03898 
0.43216 
0.48873 
0.04013 


P7,7 


9 
800 


0.01125 


0.01178 


P7,5 


39 
160 


0.24375 


0.24244 


P7,3 


463 
800 


0.57875 


0.57933 


P7,l 


133 
800 


0.16625 


0.16645 



Table 1: Calculations of PN,k, the probability of finding k real eigenvalues from an N x N 
matrix Y = A^ 1 B. The second column is the analytic calculation, the third column is the 
analytic calculation in decimal. These are compared to the right column, which contains 
the results of a numerical simulation of 100,000 matrices. 
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4 Correlation functions 



We would like to make use of knowledge of the Pfaffian form of the generating function ( 29 ) 



and the skew-orthogonal polynomials (37), to compute the correlation functions P( kl ,k 2 )- 



The latter specifies the probability density for k\ eigenvalues occurring at specific points 
on the unit circle, and £?2 eigenvalues occurring at specific points in the unit disk. Note 
that there is no conditioning on the number of real eigenvalues. The probability density is 
normalised so that P(k 1 +i,k 2 )/ 1 P{k x ,k 2 ) corresponds to the density of eigenvalues at a specific 
point on the unit circle, given the location of the k\ eigenvalues on the unit circle already 
specified, and the ky eigenvalues in the disk already specified. It can be calculated in terms 



of the summed up generalised partition function ( 28 ) by functional differentiation, 



P(fc X) fc 2 )(e,w) 



1 



-k 2 



-Z N [u,v] 



u=v=l 



(53) 



Zn[u, v] 5u{e\) ■ ■ ■ 5u(ek 1 )Sv(wi) ■ ■ ■ 5v(wk 2 
To compute (53) from the formula (29) for Zjy[u,v], we draw on established theory 



relating to calculation of ptk lt k 2 ) f° r the rea ^ Ginibre ensemble. For N even, the summed 
up generalised grand partition function for the real Ginibre ensemble is proportional to 

Pf [a jik + P j>k ] jt k =i,...,N, (54) 
where, for arbitrary monic polynomials Pi(x) of degree i, 

/oo poo 
dxu(x) / dyu{y)e~^ x2+y2)/2 p j -i{x)pk-i{y)sgT).{y-x), 
-oo J —oo 

j3j )k = 2i dxdyv(x,y)e y2 ~ x2 eric(V2y) 

x (pj-i(x + iy)p k -i(x - iy) - Pk-i(x + iy)p j - 1 (x - iy)^J . (55) 

Comparing this to p9] ) and p6] ) we see that upon the identifications Pk(x + iy) qk(w), 
Pk(x — iy) o q k (l/w), the two expressions are structurally identical. In the case of (54), 



with {pk(x)} chosen to have skew-orthogonality properties analogous to (36), a (ki + /c2 



(ki + ^2) Pfaffian formula for P(k lt k 2 ) nas been deduced which makes use of the structural 
properties exhibited by (55) (see [9]); this can be adapted to the present problem, allowing 
us to deduce the following explicit evaluation. 

Theorem 4.1. Let D denote the unit disk in the complex plane, and let dD denote its 
boundary, the unit circle in the complex plane. For N even 



P(fci,fe)( e > w ) = Pf 



Kjsr(ei,ej) K N (ei,w m ) 
K N (wi,ej) K N (wi,w m ) 



ei G dD, Wi G D, 



K N (s,t) 



D(s,t) 
-S(t,s) 



S(s,t) 
I(s,t) 



(56) 



(57) 
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with 
where 



and 



D(xi , Xj ) 

S {x i , Xj ) 
Xj) 



dj (x) 
bj{x) 
e(xj , Xj ) 

n 



e = {ei, ...,e kl } , w = {wx, -, w k2 } 

2 1 r 
' — » 77 L 

z=o 

JV_1 

i r 

^ — [02/(^)^+1 (^i) - a2i+i(xi)b 2 i(xj) 



-1 



1 r 1 

^ - \b 2 i(xi)b 2 i + i(xj) - b 2 i + x(xi)b 2 i(xj) +e(x il x j ), 

1=0 ri 



\x\ 1 r(x)qj(x), x G D, 
y 7 — i/2 r(x)qj(x), x G 

|x| _1 r(x _1 ) (x -1 ), 

/ o 2u r(e i9 )g j (e i0 )sgn(0 - arg(x))d#, x G &D, 

sgn(arg(xj) - arg(x,)), Xj, Xj G <9D, 
0, otherwise, 

ai + Pi, 



and the polynomials {qi(x)} are as in (31). 



Due to the polynomials of Proposition 3.7 being skew-orthogonal in the sense of (50) 
we can adapt the results of [35] to the present problem to yield the correlations for N odd. 



Theorem 4.2. For N odd, the correlations obey the same Pfaffian form as in (56) above, 
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with the kernel structure as in (51). The kernel elements are given by 

r(jv-i)/4i-i l 

D(xi,Xj)= — a 2 i(xi)a 2 i+i(xj) - a 2 i+i(xi)a 2 i(xj) 



1=0 



n 



(iV-l)/2-l 



+ E 



i=[(AT-l)/41 



a 2 i(xi)a 2 i + x(xj) - a 2 i +1 (xi)a 2 i(xj) 



r(Ar-i)/4l-i 



(jV-l)/2-l 



+ E 



K(tf-i)/41 



Z=0 



n+1/2 



?7 



S(xi,xj)= ^2 ~ a 2l(xi)b 2 i+i(xj) - a 2 i+i(xi)b 2 i(xj) 



a 2 i(xi)b 2 i + i(xj) - a 2 i +1 (xi)b 2 i(xj) +K(x il x j ), 



r(iv-i)/4i-i 



1=0 



n 



I{xi,Xj)= - b 2i{xi)b 2 i +l {xj) -b 2 i +1 {xi)b 2 i(xj) 



(AT-l)/2-l 



+ E — 

i=[(jV-l)/41 r ' +1 /2 



M^O^+iOe.;) - ^a+i^i)^^') +e(xi,Xj) +a(xi,Xj), 



where a,b and e(xi,Xj) are as in Theorem 4-1 

K,{xi, Xj^) = 



otherwise, 



0, 



V^i h N-i{xi) - b N -i(xj)), Xi,Xj£dD, 

-^b N -i{xj), Xi£dD,Xj£D, 

^6Ar_i(xj), Xi£D,Xj£dD, 

, 0, otherwise, 

n+l/2 = «Hl/2 + A+l/2; 



and the polynomials {qi(x)} are as in Proposition 3.1 
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4.1 Kernel element evaluations 



Clearly, the correlations in (56) are completely determined by the kernel KN(s,t) of (57). 
The elements of the kernel satisfy the following relations 

r6 2 



I r , r (ei,e 2 ) 



/ 5 r , r (e,e 2 )d0 + sgn(0i-02), 

J0! 



d 

D r , r (ei,e2) = 7^-SV,r(ei,e2), 

rw 

I rjC (x,w) = / S r>c (z,w)dz, 

J X 

D rc (ei,w) = r-^S^ciei,™" 1 ), 
\w\ A 

I c ,c(w 1 ,W2) = t—^S c>c (Wi 1 ,W2), 
\W\\ 

D cc (w 1 ,w 2 ) = T—r^S^ciwi,^ 1 ), (58) 

\W2\ Z 

where the subscripts denote the real or complex nature of the two arguments. It is clear 
from these relationships that the various S(s,t) determine the nature of the kernel block 
(57). These S(s, t) can be written in a summed-up form which is independent of the parity 
of N. 



Proposition 4.3. The elements of the correlation kernel (51) S r ^ r (s,t), S TiC (s,t), S c , r {s,t) 
and S CtC (s,t), corresponding to real-real, real-complex, complex-real and complex- complex 
eigenvalue pairs respectively, can be evaluated as 



r((iv + i)/2) 
^ (ei ' e2) = 2^r(iv/2) cos 



N-l 
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Sc,r(w>,ei) 



i\ 1/2 1 iJV /r((iV + l)/2) 

^2 N ^\j r(iv/2 + i) 

' e -i(e ro -ei)/2 e i(0«,-0i)/2 \ iV_1 



2 



I 1/2 1 -1/2 

-A 1/2 1 N(N-l) 
0fj ~r~ w 2*+ 2 ^ 



(ft 



(1 + t2)7V/2+l 



1/2 



1/2 



] X -r w (1 + *2)JV/2 



+1 



,-i{6x-0 w )/2 e i(6i-e w )/2 
+ ~ 



N-2 



,-i{6 1 -e w )/2 



/ r((jv + i)/2) 
/ r(iv/2 + i) 

,i(0i-M/2\ 



r, 



1/2 



-1/2 



S c>c (w,z) 



N(N — 1) 



oo 
-1 



(1 + i 2 ) 



JV/2+1 



1/2 



,1/2 

oo 
-1 



r u 
dt 



1/2 



(!+*») 



7V/2+1 



1/2 



e i(0z-0w)/2 e -i(e z -e w )/2 
+ 



N-2 



J(9 z -e w )/2 



-i{0z-O w )/2 



\ (r w r z ) 1/2 {r w r z ) 1 / 2 J 

where w, z := r w e t8w , r z e l9z . 

Proof: Using the binomial theorem, the identity 



{rwrz) 1 / 2 (r w r z ) 1 / 2 



n-l 
p=0 



In- 1 
2p 



~2p-l 



and the results of Proposition 3.2 (for the even case) and Proposition 3.7 (for the odd case) 
the respective sums can be performed. □ 
Note that J^ n S r>r (e-\_, e\) d6 = En, providing a further derivation of (Ji). 

The simplest cases of Theorem 4.1 are (fei,^) = (1>0) and (fci,/^) = (0,1). Since 
these correspond to the real and complex densities respectively, we write P(i.o)(e) = P^(9) 
and P(o,i)( w ) = P(i)( w )- According to Theorem 4.1, p^(0) = S r>r (e,e) and pfLJw) = 
S CiC (w,w), and we read off from Proposition 4.3 the evaluations (|8|) and (17) respectively. 
Recalling Q has the large iV form 



P\i)(0) 



while integration by parts of (17) shows 

(N-l) 1 



p(i)H ~ 



1 ttN 
27V^~' 



N — 11 



7T 



(1 + r 2 ) 2 iY-2vr(l-r 2 ) 2 



+ 



N 
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valid for r € [0, 1 — 0(l/vJV)]. To leading order in N the eigenvalue density is therefore 
equal to 



N 1 

7r (1 + r 



2 , 2 (59) 



for all r £ [0,1]. This, projected stereographically onto the half sphere, gives a uniform 
distribution. This 1/N convergence should be contrasted with the exponential convergence 
in the case of the polynomials ^ (see |31j). 

4.2 Averages over characteristic polynomials 

As emphasised in |19} |3| there is a large class of eigenvalue jpdfs such that the eigenvalue 
density is given in terms of an average over the corresponding characteristic polynomials. 
This is true of the one-point function (density) for the complex eigenvalues, with N — > N+2 
(for convenience), in the present generalised eigenvalue problem for which the jpdf is given 
by ([9}. Thus write 

k (N+k)/2 

C N (z) = l{(z-e j ) H (z-w s )(z-l/w s ) (60) 

j=l s=k+l 

for the characteristic polynomial in the N x N case ofY = A^ 1 B conditioned to have k 
real eigenvalues, with eigenvalues transformed according to ^ and ([7]). Letting 

Gn := t ^S^ n(N + l)/2)^r(iV/2 + 1)^ f[ 



2 N(N-l)/2 i II") ll r (s/2) 2 ' 

s=l 



which is the pre-factor in (10), then it follows from ^ and the definition of the density 
that 



p (N+2,c) (z) = 1 T{z)T{1/m /- z _ z){Cn (z)C n (1/z)), (61) 

G N \z\ l 



where the superscript iV + 2 denotes the number of eigenvalues in the system. Of course we 



can therefore read off from ( |17| the exact form of the average in (61 ). Moreover, in keeping 
with the development in [3], we can use our integration methods to compute the more 
general average (Cn(zi)Cn(1/z2)) which we expect to be closely related to S c>c (zi, Z2), in 
accordance with known results from the real, complex and real quaternion Ginibre ensemble 

123 HID El M- 



Note that we will also introduce a superscript on a and j3 to indicate the size of 
tern that the 
normalisations. 



system that they relate to, that is a^j has j, I = 1, t, and are the corresponding 
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Proposition 4.4. With the characteristic polynomial Cn given by (60) and (■) an average 
with respect to |5p summed over k, one has 



(C N ( Zl )C N {l/z 2 )) 



Gn 

Gn+2 
N/2 



(l/*2 " Zl) 



1 



-1 



=0 a s + pi 
Gn 



Gn-\ 



where it is assumed N is even. Furthermore 

G N +2 t{zx) t(1/z 2 ) 



V \ z l\ \ z 2\ J 



N^N+2 



(62) 



S c ,c(zi, z 2 )\ 



N^N+2 



G N \z\\ \z 2 \ 



(l/^-«i)(Cj>r^i)CW(l/f2)) 



from which we reclaim (61). 
Proof: From ^ we see that 



(N+k)/2 



v-1 







' 1 > 


Ws 






1 






w 


Zl, 


1) 



Integrating over e and w gives 

(C N ( Zl )C N (l/z 2 )) = G N {l/z 2 - zx)' 1 

/c fixed 

x [x^ 2 ][ C fc ]Pf [x(C 2 a^+ 2 > + ^) + ^ 

where 7]*/ = qj-i{zi)qi-x(l/z 2 ) - qi^i(z 1 )q j - 1 (l/z 2 ) = l,...,t), and the parameters 
u, v in a[u, v] and f3[u, v] are taken to be one. Summing over k leads to 

N 



(C N ( Zl )C N (l/z 2 )) := (C N (zi)C N (l/z 2 )) 



k=0 
k even 



k fixed 



G N (i/- Z2 - zl y^ N/2 n x(*^ +2) + 4r 2) ) + t17 +2) 
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Using the skew-orthogonal polynomials (37), we find 



^ 2 ]Pf [>c{*f l +2) +^ +2) )+^ +2) 



N/2 N/2 
^ 72s+l,2s+2 \\{ 



(N+2) 3 (N+2) , 

"2j+l,2j+2 ^ ^2j+l,2j+2J- 



3=0 



The evaluation now follows upon recalling the form of Zjv+2(1) = 1 = Gn+2 TiiHo ( a i + 
p(N+2)^ j m p^ e( ^ D y p r0 position 



3.3 



The expression for S cc (zi,Z2) 



AT-s-iV+2 



is a simple manipulation of (62). 



□ 



4.3 Scaled limit 

Before implementing the fractional linear transformations ([6]) and ([7]), we have from |3l) 
that the density of real eigenvalues near the origin is proportional to En, and thus y/N. 
A scaled limit involves changing the variables so that this density, and that of the com- 
plex eigenvalues, becomes of order unity. Such a limiting procedure is of interest because 
the resulting correlations are expected to be the same as for the generalised eigenvalue 
problem with entries chosen from general zero mean and finite variance distributions, and 
furthermore the same as for the eigenvalues of the real Ginibre ensemble, scaled near the 
origin. 

In the complex case, an analogy between the eigenvalue jpdf of the generalised eigen- 
value problem and the Boltzmann factor for the two-dimensional one-component plasma 
on a sphere [10J, together with the analogy between the eigenvalue jpdf for the Ginibre 
matrices and the two-dimensional one-component plasma in the plane \i>, allow this latter 
point to be anticipated from a Coulomb gas perspective. 

The limiting correlation of the eigenvalues in the vicinity of the origin for the real 
Ginibre ensemble have recently been computed in [9j Corollary 9] (see also p21 [36] ) . 

Proposition 4.5. For random N x N real Ginibre matrices the correlations for the eigen- 
values in the vicinity of the origin are, in the limit N — > oo, given by 



Pfei,fc 2 (x,z) = Pf 



[K r)r (xj,Xi)} jt i = i^. M [K r>c (xj,zi)] 



j=i,...,*i 



-K r , c (xj, zi)] 

i=i,...,fca 



T 



[K c , c {zj, Zf)]i,i=i,-,*2 



(63) 
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with 



K r>r (x,y) 

K rtC (x, w) 



i_ e -(x-y) 2 /2 



erfc(\/2Im(2;)) 



J_ P -(*-2/) 2 /2 



2vr 



hsgn(x — y)erfc 



\ x ~y\ 

V2 



(w - x)e~^~ w ^/ 2 i(w - x)e- { - x -^ 2 / 2 

_ p -{x-w) 2 /2 _j p -{x-w) 2 /2 



K C)C (w, z) 



0,7V 



erf c ( V^Im (w ) ) er fc ( V2Im ( z ) ) 



-(w-z) 2 /2 



w e 



2; 



- w)e-^-~ z ^/ 2 



In the present problem, with our use of the transformed variables ^ and (J7]), the 
original origin has been mapped to (1,0). We must choose scaled co-ordinates so that in 
the vicinity of this point the real and complex eigenvalues have a density of order unity. 
For the real eigenvalues, from the knowledge that their expected value is of order \/N and 
that they are uniform on the unit circle, with e,- := e lXj , we scale 



2X 



Xj H-> 



N 



(64) 



For the complex eigenvalues, which total of order N in the unit disk, an order one density 
will result by writing 

Wj 1 + ^LW,. (65) 
3 s/N 3 

Note that the real and imaginary parts have been interchanged to match the geometry of 
the problem in the Ginibre ensemble, that is so the eigenvalues are again distributed in the 



upper half-plane, including the real line. The factors of 2 in (64) and (65) are included so 



an exact correspondence with the results of Proposition 4.5 can be obtained. 

(x, x) is interpreted as a density, the normalised quantity is S r>r (x, x)dx. It fol- 
lows then that in the more general case we must look at the scaled limit of S rr (x, y)\/dxdy 
and S CtC (wi,W2)Vd 2 w\d 2 W2- For S rfi (x,w) and S C:r (w,x) we require that the product 
S rjC (x,w)S C;r (w,x)dxd 2 w has a well defined limit. From (64) and (65) we see 

V ax ay 1— > -^VdXdY, 

Vn 

1 4\ 3/2 
dxd 2 w i—)- I — J dXd 2 W, 



^d 2 Wl d 2 w 2 1 — y —Jd 2 Wid 2 W 2 . 



N 
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With this change of variables the large N form of the correlation kernel for the spherical 
ensemble matches that of the Ginibre ensemble. 



Proposition 4.6. Recall Kjy(s,t) from (51). Replacing Xj and Wj according to (64) and 
(65) then taking N — > 00 gives 

1 

K N (ei,ej) ~ K rr (Xi,Xj), 



'N 

2 5/2 

V2 



K N (ei,Wj] 
K N (wi,ej] 



N 



K N (wi,Wj 



-(K r , c (Wi,Xj)) T 
K c , c (W h W 3 ). 



Proof: From the explicit functional forms of Proposition 4.3, we see that elementary 



limits suffice. For example, changing variables 1 1— > 2t/y/N shows 

dt 



(1 + t2)7V/2+l y 2N 



7T 



erfc(\/2ImTy). 



Combining such calculations we obtain 



'N 

2 5/2 



1 



N 
V2 



2tt 

i 



S C)r (wi,ej) ~ 



2?r 

i 



erfc(v / 2ImW J ) e~ {Xl ~ w ^ /2 i(W j - X { 



'2-k 

1 



erfc(^/2ImWi) » e -( w *-- x 'i) 3 /2 ) 



erfc(\/2ImWi) J ertc(y/2JmW s 



xiiWj-W^e-^-^ 2 , 



which is in agreement with the off-diagonal entries on the RHS of the present proposition, 
as implied by Proposition 4.5 (when one recalls that <SV, C and S c , r never appear individually; 
only as the product 5 rjC 5 c>r ). 

Recalling the inter-relationships (58), the other kernel elements D and I can be obtained 
from S, giving the diagonal entries required by Proposition 4.5. □ 
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4.4 Sum rules 



With A as in ([9]) we have 
A(e,w,- 



cxp 



k (N+k)/2 

Yl l °g\ e p- e j\-Yl Yl log 

i<i<p<fc i =1 s=k+i 



1 1 


M) 


\w b 


W a \) 



Yj \o&\w h -w a 

k+l<a<b<(N+k)/2 

This is the Boltzmann factor of a two-component log-potential Coulomb gas, consisting 
of k unit charges at {&j}j=i,...k confined to the unit circle, (N — k)/2 unit charges at 
{t«o}o=fc+i,...,(JV+fc)/2 confined to the unit disk, and a further (N — k)/2 image charges 
to those in the unit disk, which are at positions {l/wa}a=k+l,...,(N+k)/2 outside the unit 
disk. The other factors in ([£]) are, from a Coulomb gas perspective, one-body terms due 
to the coupling of the charges to an external background charge density. We have seen 



in Proposition 4.6 that in a certain scaled limit the correlation functions for this two- 
component Coulomb gas tend to functional forms known from the study of the real Ginibre 
ensemble. These should exhibit features characteristic of a two-component Coulomb gas. 
Here we will exhibit two such features of the scaled correlations. 

In regard to the first of these, suppose that in the limiting system we fix real eigenvalues 
at {xj}j = i v .. ) fc 1 and complex eigenvalues at {zj}j=i,...,k2 (the latter also requires complex 
eigenvalues at {%}j=i,...,fc 2 )- Regarding this action as perturbations, the Coulomb gas 
perspective tells us that to maintain equilibrium the system will respond by surrounding 
the fixed eigenvalues with a screening cloud equal and opposite in total charge to that of 
the perturbation. In terms of the correlations this gives rise to the sum rule (see e.g. |144 
eq. (14.20)]) 



/oo 
pffei+i,^)^^}^!,-^! u ivh i z j}j=i,..,k 2 ) d v 
-oo 



+2 



pfk u k 2 +i)({ x j}j=i,-,k^ {zj}j=i,...M U {z}) d 2 z 



-(ki + 2k 2 )pf klM) ({x j } j=h ... M ;{z j } j=h ... ik2 ), 



where z = X + iY, d 2 z = dXdY, and pj£ k2 }({xj}j 



(66) 

i,...,fci) { z j}j=i,...,k 2 ) ^ S *he truncated 
correlation function (see |14|, Eq. (5.3)]). The factor of 2 with k<i in (66) is due to the 
complex eigenvalues always occurring in complex conjugate pairs. 

The explicit form of the pT ki k ^ can be read off from (63). For this it is convenient to 
introduce the quaternion determinant, qdet, according to 



qdet AZ 2N = Pi A, 



Z 2 n '■= ljv ® 



-1 

1 
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for A a 2N x 2N antisymmetric matrix. The correlation function (|63|) then reads 



[K rr (Xj,Xt)]ji = i fa [Kr,c( x ji z l)]i =1 —> k i 

! = l,...,fc 2 

[K C} r(zi,Xj)] i=l,...,fc 2 [K CiC (zj, Zl)]jJ = i^k 2 



j=l,...M 



where, with K c ^ r (z,x) := — (K r:C (x,z)) T , each 2x2 block K*(u,v) is related to the corre- 
sponding block in (|63|) by 



K*(u, v) = K*(u, v) 



-1 

1 



The advantage of such a determinant form is that it then follows (see e.g. [14, Eq. (7.184)]) 
that the corresponding truncated correlations are given by a sum over maximum length 
cycles in the determinant, 

PfktM) ( x > z ) = (-l)^ 2 " 1 (K( yil ,y t2 )K(y i2 ,y i3 )...K( mki+k2 , yil )y 0) . (67) 

cycles 
length 

Here {yi}i=i,...,k 1 +k 2 = {»i}i=i,...,jfc 1 U {zj}j=x r ..,k 2 , and the operation (.)(°) refers to ^Tr. 

The sum rule (66) is a corollary of integration formulas involving the product of two 
matrix kernels K*. 

Proposition 4.7. We have 

K r r (x,u)K r r (u,y) du + 2 / K r c (x, z)K c r (z, y) d 2 z 
> Jr 2 + 

T > , x [o 0l [1 ol ~ . , 

Kr,r{x,y) „ , + „ n K rtr (X,y), 



k c>r {z,y)K r>c (y,v)dy + 2 / KcAz, w)K C:C (w, v) d 2 w = 2K c>c (z, v), 

Jr 2 + 

K cr (z,y)K rr (y,x)dy + 2 / K CtC {z,w)K c , r (w,x) d 2 w = K CtT (z,x), 
K r ^ r {x,y)K rfi (y,z)dy + 2 / K r , c (x,w)K CiC (w,z)d 2 w = K r ^(x,z). 



Proof: Each of the above equations requires evaluating the integrals for the four entries 
of the matrix products. We will illustrate the required working by giving the details in the 
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case of the (ll)-component of the first of the equations. With the notation (X)jk denoting 
the entry (jk) of the matrix X, we have 

K r , r (x,u)K r , r (u,y)du) n = —J e-^-^-^-y?du 

^= / (x - u)sgn(u - y)e~^ x ~ u ^eric( — 
'2tt J-oo ^ v2 ' 



2 

Completing the square shows 



du. 



j_ r i (a! -„)2_i( U _ V )2 du = j_ -!(*_„)» 

2^ Loo 

while writing (x — t^e - ^"") = Jj-e - ^ - ") 5 integrating by parts and making further use 
of the above integral evaluation shows 

— \= f (x - u)sgn(u - y)e^^ x ~ u ^eifc ( ^— — du 



1 =e -2(^) 2 + — e -\{x-v) 



2 



\/2tt 2^/tt 

Hence 



■DC 



K rr (x,u)K rr (u,y)du) = -Leri^-f) 2 + ^ e -|(^) 2 . (68) 



For the integral over in the first of the equations, we have 

aK T)C {x, z)K cr (z, y) d 2 z ) 
,2 / 11 

= Re- f ^c{^j2Y){z - x)e-^ x -~ z)2 -^ y - z)2 dXdY. 
k J«. 2 + 

Recalling that z = X + iY, completing the square in X and Y, translating the integral in 
X, and noting Xe~ x dX = reduces this to 

Re- f erfc(V2Y)(Y - i ^y)e- x2 ^ Y -^ 2 dXdY 



Now writing \ Y — i x ^-je^ Y 1 2 V ) = ±JL e ( Y 1 2 V ) anc l integrating by parts allows this 
integral to be evaluated, and we obtain 

K r , c (x, z)K c , r (z, y) d 2 z) = - J_ e -(^) 2 /4 + _L e -^-v?l\ ( 69 ) 
/ 11 ^yvr \/2vr 
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Forming (|68|) plus twice (69) gives ^J^e ( a J 7 ) 2 / 2 , which is the (ll)-component of K TjT {x, y), 

□ 



in keeping with first of the equations of the proposition. 



To prove the sum rule (66) we substitute (67) for the truncated correlations on the 



LHS of (66). We see that the required integrations can be computed using Proposition 



4.7 



When involving the first of the integration formulas therein, we can either use the 
cyclic property of the trace, or sum together pairs of terms, to effectively replace the RHS 
by K r ^ r {x,y). In all cases this allows the expression resulting from the integration to be 



identified with a term in the expression for —p^ ^ on the RHS of (66) as implied by (67). 
Moreover, each term is repeated (k± + 2^) times, thus verifying (66). 



In addition to the sum rule (66) there is a second sum rule satisfied by the correlation 
functions, as suggested by the Coulomb gas analogy. For this we consider the complex 
moments of the screening cloud due to a fixed complex eigenvalue at point z. This screening 
cloud is defined as the function of w £ R 2 . and x £ K given by 



2 / ? fo,2) 0, w )x w m 2 , + pji,i) i z > x )Xxm + 25 (2) (z - w)p (0jl) (w). 



(70) 



We know as a special case of (66) that integrating this over w and x gives zero. Another 
feature of the screening cloud, seen by inspection of Proposition 4.5 is that for large \w\ 



and \x\ it decays at a Gaussian rate to 0. In the theory of Coulomb systems (see e.g. |30j) 
this rapid decay can be shown to occur only if the complex (multi-pole) moments of the 
screening cloud all vanish. We know from |24t IT5] that this statement must be modified 
when image charges are present so as to relate to the screening cloud of the charge/image 
system. Consequently, we should re-interpret (70) as the function of w £ M. 2 \R. and i£l 
specified by 



P{o,2) (*> w )x we R 2 + + pfi,i) ( z > x)xxm + <5 (2) (z - «;)p(o,i) (H 



(71) 



with the property that 



P(o,2)( z > w ) = £(0,2)0;™) P(o,i)(w) = P(o,i)H- (72) 
By translation invariance of the system in the x-direction we can set z = iyo- We then 



see that with the conditions ( 72 ) the odd complex moments of ( 71 ) vanish by symmetry. 



whereas the vanishing of the even complex moments requires 



w 



2p Pio,2)( z ^ w ) d2w + I x 2p pf 11) (z,x)dx = -2z 2p p {0:1) (z) p = 0,l, 



Recalling the explicit form of the truncated correlations as implied by ( 67 ) and Proposition 



4.5, by multiplying both sides by a p /p\ and summing over p we obtain an equivalent form 



of this sum rule, which we state and prove in the following result. 
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Proposition 4.8. Let z = iyo and suppose \a\ < 1. We have 
2 [ e aw ' 2 K C)C (z, w)K C;C (w, z) d 2 w + f°° e ax ' 2 K c>r (z, x)K r , c (x, z) dx = 2e az2 R CjC (z, z). 

JR 2 + J-oo 

Proof: We will illustrate our methods by considering the (ll)-component. For the first 
term on the LHS this component can be written 



-erfc (V2y ) / dY dX e a{x+iY) erfc(v^y) 
t7r Jo J-oo 



Separating the terms involving the partial derivatives, and integrating by parts in each 
reduces the double integral to 

/oo 
e aX2 ~ x2 dX 

8 ^ dY j°° dX (ia(X + iY) - Y y^ x +^) 2 ^ Y2 ~ x2 {e 



a(X+iY) 2 -2Y 2 ~X 2 (p (Y+y y _ e (Y-y ) 2 



Upon completing the square in X, simplifying, then completing the square in Y the second 
of these integrals can be evaluated, giving that the (ll)-component of the first term on the 
LHS is equal to 

1 . (•OO . 

—erfc (V2y ) ( - 4y e y ° J e -( 1 ~ a ) x2 dX + 8V&re 2y2 - ay2 y ) . 

An analogous strategy in relation to the (ll)-component of the second term on the LHS 
shows that it is equal to 

oo 



1 A 00 

-erfc (V2y )y e y2 / e~ {1 ~ a)x2 dX 
n J-oo 



Adding together the above evaluations of the (ll)-components of the terms on the LHS 
gives 



- ay o{2^evtc(V2y )y e 2y o), 



which we recognize as the (ll)-component of 2e ay o K c ^ c (z, z), as required by the RHS of 
the sum rule. □ 
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Appendix A 



The following lemmata from [33] are required in the proof of Proposition 2.1 



Lemma I (Theorem 2.1.5). For X = aY{3, where ap x p and (3qxQ are arbitrary real 
matrices and Yp X Q has PQ independent entries (ie. the wedge product (dY) has PQ 
factors) then 

(dX) = \det(a (3 T )\ (dY) 

= |det(a) Q det(/3) p | (dY). 

Lemma II (Theorem 2.1.14). For any N x M (M > N) matrix X, ifW = XX T then 

(dX) = c detW (N ~ M - 1)/2 (dW), (73) 

where c is independent of W. 

Lemma III (Theorem 2.1.6). For A an N x N real non-singular matrix and M , an N x N 
real symmetric matrix, one has 

(A T dMA) = det(A T ^)( JV - 1 )/ 2+1 ( ( iM). 

Lemma IV (Ch.3, Equation (22)). For C a real symmetric N x N matrix 

N 

(dC) = ~Y\_( x k ~ Xj)dx\ ■ ■ ■ dxN(R T dR), 

j<k 

where the {xi} are the ordered eigenvalues of C , and R is the real orthogonal matrix of 
eigenvectors. 



A.l Proof of Proposition 2.1 



Recall that Y = A l B. By letting a = A and /3 = ljv in Lemma ||] we see that 

(dB) = \detA\ N (dY), 



so we rewrite (18) as 



V(Y)(dA)(dB) = (27T)- N2 e- 1 2 T < AAT ^ +YY ^)\detAA T \ N / 2 (dA)(dY). 

(Here and below (dA)(dB), and similar, is to be interpreted as the wedge product of 
the corresponding differentials.) Setting C := AA T , Lemma [n] tells us that (dA) = 
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c(detC) 1 / 2 (dC). Integrating over C (noting that C is positive definite, denoted C > 0) we 
have 

V(Y)(dY) = (2ir)- N2 c [ (detC^ N -^ 2 e- k M^+yy T ))(dC)(dY) 

Joo 

= {2k)- n \ [ (detC)( iV - 1 )/ 2 e -5 Tr ((^+^ T ) 1/2 ^^+^ T ) 1/2 )(dC)(dy). 
Joo 

Carrying out the change of variables C -> {1 N + YY T ) 1 / 2 C(1 N + YY T ) 1 / 2 we use Lemma 
EH to find 

V{Y){dY) = (27r)- N2 cdet(l N + YY T r N [ det^^-^e-^^^C)^). 

Joo 

Using Lemma [IT] we can calculate c according to 

J e -Tr(AAT)/2 {dA) = c f e -Tr(C)/2_ (dC) 

and so 



Joo det(C)V2' 
(2^)^/2 



telling us that 



/ c>0 det(C)-V2 e -TKC)/2 ((ic) - 



= (2,)- ^det(i„ + ry' )- ^ d< : t(c) - 1/2e ^, c)/2(dc) (dY). 



bmce C = AA T is symmetric, using Lemma |IV| the ratio of integrals can be rewritten as 

J c>0 det(C)( Ar - 1 )/ 2 e-i Tr ( c )(dC) 
f c>0 (detCy^e-^C)/2 {dC) 

_ J(o,oo)iv IlHi ^ (jV ~ 1)/2 e- x '/ 2 nf< fc - Sj|<fai — dxN 

/(o,oo)w n£=i ^r 1/2e_a ' i/2 nf< fc kfc - ^i^i ■ • • rfxjv 

which is seen to be a ratio of Selberg-type integrals, which have known evaluations in terms 
of gamma functions (see e.g. [HI Ch. 4]) . The result now follows. □ 

Appendix B 

The purpose of this appendix is to integrate over Rn, the strictly upper triangular elements 



of i?jv in (19), leaving us with just the dependence on the eigenvalues. This will be done 
column-by-column, starting with the case where j > k, that is, the columns corresponding 
to the complex eigenvalues, and then proceeding onto those columns corresponding to the 
real eigenvalues. 
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B.l Complex eigenvalue columns 

In the region j > k of Rn we can isolate the last two rows and columns to write 



Rn 



Rn-2 U 



T 



where u is of size (N — 2) x 2 and T is of size 2 x (N — 2). So then 



Ijv + RnRn 
and 



l N -2 + Rn-2R N ^ 2 + UU ~ UZ^fa + ZmZ^) Z m U 



1 2 + Z m zl 



det(l7v + RnR T n) 



det(l 2 + ) det(ljv_2 + Rn-2Rn-2) 

x det(l 2 + (1 2 + 2; m ^)~ 1 n T (lAr_ 2 + Rn-2Rn-2)~ 1u ) 



where we have used the general identity (for appropriate sized 1) 

det(l + AB) = det(l + B^). (74) 
We are now in a position to integrate over the elements of the matrix u 

r {du) _ i 

J det{l N + R N R%) N ~ dct(l 2 + z m zl) N det(ljv_ 2 + R N ^ 2 R T N _ 2 ) N 

x r m 

J det(l a + (1 2 + z m z^)-^u T {l N ^ + R N - 2 R T N _ 2 )^u{l 2 + z m ^)-i/2)iv . 

where the integral for each independent real component of u is over the real line. Changing 
variables v = (ljv-2 + Rn~2RJj^2)~ 1 ^ 2u (^-^ + z m z m) -1 ^ 2 we use Lemma [i to find 

(du) 1 



det(l N + R N R%) N det(l 2 + z mZ T)m+i det(l A r- 2 + R N - 2 R T N _ 2 ) N ~ l 

(dv) 

2 + v T v) N ' 



I det(l 



Iterating over all columns corresponding to complex eigenvalues we have 
(du k+1 ) 1 



r (du N - 2 ) ■ ■ ■ 

J det(l7v + i? 



n 



R N R T N ) N det(l k + R k R T k )(^)/z 

I f (dv N - 2 -2s) 



n 



det(l 2 + tf «- - - 



(75) 
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where the subscripts * on the matrices v*,du*,dv* denote their number of rows. 

To evaluate each of the (AT — k)/2 integrals we use a similar method to that used in 
Proposition 2.1 Firstly, for each vn-2-2s, we let vJ l _ 2 _ 2 s v N-2-2s = C and apply Lemma 
|II]to get 

and 



(do) = c(detC) (JV - 2s - 5)/2 (dC) 



(76) 



cj{&etC)( N - 2s -^ 2 e- T < c \dC) 

And so, with k := (A - 2s - 5)/2, 
(dv N - 2 -2s) 



N-2s-2 



JV- 



det(l 2 + %_ 2 _ 2s f at-2- 



7T 



7T 



iV-2s-2 



J N-2-2s l 
roo poo 



-2s) 



X- 



JO 

oo roo 



(l + ar 2 ) 



-25-2 / (det C) K det(l 2 + C)*- N jdC) 
/(det C) K e- TrC (dC) 



'0 «/ 

_AT-2s-2 



7T 



cc'fa^e " ; e ~*|a;i — xi\dx\d%2 

1 /•! 



J 

oo coo 



JO 



vxy K 2^ - yi) (JV ~ 1)/2 (i - y 2 ) (w - 1)/2 |yi - wl^i/a 

-1 



x^x 2 e 11 e 2:2 1 xi - x 2 |dxi<ixi 



(77) 



where use was made of Lemma IV for the second equality, and the change of variables 
y = x/(l + x) for the third. 



find 



A ratio of Selberg integrals has again appeared in (77), and using results from [T3] we 

v-2,-2 r((iv + i)/2) r(Jv/2 + i) 



J det(l 2 + v T N _ 2 ^ 2s v N ^ 2s )- {N ~ s \dv) = tt 1 



r(JV-a-l/2) T(N-s) 



(78) 



The case A odd, k = 1, corresponding to s = (A — l)/2 — 1 is special since then vn-2-2s 
consists of 1 row and 2 columns, and thus is the only case in which the number of rows is 
less than the number of columns. We must then write 

det(l 2 + V%_ 2 _ 2s V N - 2 -2s)~ P = (1 + V N -2-2svJ]-2-2s)~ P i 



using (74). However, it turns out that the change this implies to (77) does not effect the 
evaluation (78), even though (76) is no longer valid. So in all cases, after having integrated 
over the Rij for % > j in the columns corresponding to complex eigenvalues we are left 
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with 



(duN-2) ■ ■ ■ (du k+ i) _ (Ar ^. )/2 ] 



det(l N + R N R%) N det(l 2 + z s zJ) N /i+i 

-*)/ 

n 



, '' Y " / ''' ^\n-2 S -2 r((JV + i)/2) rCjv/2 + i) 1 

r(JV-a-l/2) T(iV-a) det(l Jfe + 22 fc l#)( iV +*)/ 2 ' 



s=0 



It remains to compute the integrals over the columns corresponding to the real eigen- 
values. 



B.2 Real eigenvalue columns 

We see that we are left with a function of R k , which is the upper-left sub-block of Rn and 
we isolate the last row and column 



Rk 



Rk-l Uk-1 



where now u k -\ is of size (k — 1) X 1 and T is of size 1 x (k — 1). Following the same 
procedure as for the complex eigenvalue columns, we find 

det(l fc + R k Rl) = (1 + X\) det(l fc _i + i^-i^Li) 

x(l + (1 + Agr^iClfc-i + R k - X R£_ x )- X Uk-\). 

Setting Vk-i = (lj-i + i?j-ii?J_ 1 ) _1 ^ 2 Ufc-i(l + A 2 ) -1 / 2 and again making use of Lemma|l] 
we have 



/ 



(du k . 



det(l fc + R k Rl)( N + k )/ 2 (1 + A 2 )^ 1 )/ 2 det(l fc _! + R k -xR^ N+k -^ 2 

(dv k -i) 



(l+^_l^-i) (JV+fe)/2 ' 
Iterating over the remaining columns of gives 

(eZttfe—i) • • • (dui) -A- 1 TT / 



det(l fc + R k RT)( N + k V 2 n (i + A 2 )(iv+i)/2 n y (1 + t£_ -UA _,)(2V+fc)/2-(«-i)/ 



(cf. (75)). To evaluate the integrals, we use the same method as for the integrals in (75) 
(involving Lemma |n] and now one-dimensional case of the Selberg integral, which is the 
beta integral). This gives 

(dv k -s) ,k-s)/2 r((JV + l)/2) 



7T 



(1 + v t Vk _ s )(N+k)/2-( s -i)/2 r((JV + fc - a + l)/2) ' 
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and so 



(r/ ( „_L)---u/»i) rr i r((iV + i)/2) 



fl/u \2\(JV+l)/2 II 



J det(l k + R k Rl)( N + k )/i ^(l + A^Wli r((iV + fc- S + l)/2)" 
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